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Abstract 

The Lorentz gas is a billiard model involving a point particle diffusing deterministically in a periodic 
array of convex scatterers. In the two dimensional finite horizon case, in which all trajectories involve col- 
lisions with the scatterers, displacements scaled by the usual diffusive factor \/t are normally distributed, 
as shown by Bunimovich and Sinai in 1981. In the infinite horizon case, motion is superdiffusive, however 
the normal distribution is recovered when scaling by Int, with an explicit formula for its variance. 
Here we explore the infinite horizon case in arbitrary dimensions, giving explicit formulas for the mean 
square displacement, arguing that it differs from the variance of the limiting distribution, making con- 
nections with the Riemann Hypothesis in the small scatterer limit, and providing evidence for a critical 
dimension d = 6 beyond which correlation decay exhibits fractional powers. The results are conditional 
on a number of conjectures, and are corroborated by numerical simulations in up to ten dimensions. 

1 Introduction 

The Lorentz gas is a billiard model [13], in which a point particle moves in straight lines with unit speed 
except for specular collisions with hard scatterers, most often disks or balls. It is an "extended" billiard 
in that it moves in an unbounded domain, in contrast to the usual billiard setting. Lorentz developed an 
approach along these lines to study electrical conduction in 1905 [33], however periodic arrays of scatterers 
in the plane or in higher dimensions appeared as early as 1873 (the Gallon board, or "quincunx" |21j). They 
have been a major example of hyperbolic dynamics with singularities starting from the work of Sinai in 
1970 [H]. 

Lorentz gas models are considered in the physics literature, for example arising as the molecular dy- 
namics model with two atoms and a hard potential after reduction of the centre of mass motion |17j . It is 
thus probably the simplest model of deterministic diffusion with physical attributes such as a Hamiltonian 
structure, continuous time dynamics and time reversal invariance; simpler but less physical models include 
one dimensional maps on / : K — >■ K with a periodicity property /(x -|- 1) = /(x) -|- 1 [201 and multi- 
baker maps on [0, 1]^ x Z [55]. In the former case, with / a piecewise linear map, the diffusion constant is 
a non-smooth function of parameters, as quantified in [26| . It would be interesting to extend such results 
to billiards; numerical studies of Lorentz (-like) models have suggested an irregular but somewhat smoother 
behaviour than the piecewise linear map [24l IJU [27l EH] . 

There is a substantial literature on diffusion in periodic Lorentz gases in dimension d = 2 in the setting 
of dispersing billiards, for example sufficiently smooth (C"^), convex with curvature bounded away from 
zero, and an upper bound on the free flight length between scatterers ('finite horizon'). It is known that 
the displacements of particles initially distributed uniformly (outside the scatterers) converges weakly to a 
Wiener process (Brownian motion) with variance proportional to the (continuous) time t as t — )■ cx) [7] ; proof 
of recurrence in the full space followed later [TS] HO] ■ The strategy is usually to relate the deterministic model 
to a Markovian random walk process for which the results are known or easier to prove. 

If we relax the finite horizon condition, the distribution remains Gaussian, but with a i In t scaling. The 
two dimensional case was considered by Zacherl et al. |46j who gave approximate formulas for the velocity 
correlation functions, and hence diffusion properties. Bleher [5|, proved some of the necessary results, deriving 
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some exact coefficients. Detailed numerical results for velocity correlations were obtained by Matsuoka and 
Martin [35]. Szasz and Varju [IS], completed the proofs in discrete time (n), obtaining a local limit law to 
a normal distribution, and also recurrence and ergodicity in the full space. Techniques there were based on 
work of Balint and Gouezel [Ij for the stadium billiard, which also has long flights without a collision on a 
curved boundary, and also a normal distribution with a nonstandard n In n limit law. Finally, Chernov and 
Dolgopyat [llj proved weak convergence to a Wiener distribution in the continuous time case. 

In general it is easier to prove statements about correlation decay and diffusion for the discrete time 
dynamics first; some diffusion results can then be transferred to continuous time, as discussed in Ref. [5] ITTj 
and Section [o] below. Correlation decay for billiards in continuous time have been considered [S] El |33j but 
are in general more difficult. Here we use continuous time {i) with connections to the discrete time diffusion 
process discussed in Sec. [9| 

The main purpose of this paper is to explore the infinite horizon case in dimensions d > 2, or more 
precisely, several varieties of infinite horizon that appear in higher dimensions. Recurrence is of course not 
to be expected, even for finite horizon, however a (super)-diffusion coefficient and further refinements may 
still be considered. One case in which the diffusion in higher dimensional Lorentz gases has been treated 
very recently is the Boltzmann-Grad limit [SUIISI]- Here the radius and density of scatterers scale so that 
the mean free path is constant. Here we will also consider the limit of small scatterers, but only after the 
infinite time limit. We will nevertheless be able to make connection with an asymptotic formula in [31) . 
In addition, we will find connections with the Riemann Hypothesis, along the lines of those presented for 
the escape problem of a billiard inside a circle with one or two holes [6]; this is also new for d = 2. Last 
but not least, there are some outstanding issues in d = 2 regarding the diffusion coefficients defined using 
moments and using a limiting normal distribution, which we need to clarify first (Sec. [2]); these may also have 
wider implications for our understanding of other anomalous diffusion processes. Main results are stated and 
discussed in Sec. [3| which concludes with a description of the following more detailed sections. 



2 Diffusion coefficients 

Different conventions are used in physical and mathematical literature for diffusion, and there are further 
subtleties involved with the definition of the logarithmic superdifFusion coefficient. Pick's law, dating back to 
1855, states that the flux of particles in a diffusing species is proportional to the concentration gradient. If, 
as in the Lorentz gas, the diffusing particles do not interact with each other, and the medium through which 
they diffuse (periodic collection of scatterers) is homogeneous at large scales but not necessarily isotropic, 
from physical considerations we expect a macroscopic equation of the form 



dp ^ d d 
dt dxi dx 



-D.^l^l^P (1) 



where p(x, t) is the density of particles, i and j are spatial indices, we are using the Einstein summation 
convention (summing repeated indices from 1 to d), and we denote quantities in the full (unbounded) 
space with tildes. The Dij are components of the diffusivity matrix D, which is positive semi-definite 
and symmetric. For now we assume that it is actually positive definite, corresponding to the possibility 
of diffusion in all directions. In simple cases it is often isotropic, Dij = DSij where D is "the" diffusion 
coefficient or diffusivity. The Green's function of Eq. ([I]) is a Gaussian 



'"(^'*) = ^TT!fTl^ (2) 
V47r< det D 

involving the inverse matrix D~^. The second moment of position with respect to this density is 

J p(x,t)xiXjdx — 2Dijt (3) 

The macroscopic picture may be connected to the microscopic dynamics in terms of the displacement of 
a single particle A(i) = x(i) — x(0). If we assume the microscopic dynamics is periodic with respect to a 
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lattice of finite covolume, we can define a configuration variable x modulo lattice translations, and consider 
the dynamics as a random process defined using a probability measure n on its initial position and velocity 
(x(0),v(0)). Expectations with respect to this measure will be denoted by angle brackets (). Depending 
on the details of the dynamics, it may be possible to prove that in the t — > oo limit A has the appropriate 
second moment 

and/or that it converges weakly to a normal distribution 

^^AA(0,E,,) (5) 

with zero mean and covariance matrix = 2Dij. Convergence of the process to Brownian motion (ie 
multiple time distributions) with the same covariance matrix has also been shown in some cases. Note in 
particular the factors of 2. We may want an overall diffusivity, generalising the isotropic case: 

itrD= hm -^(A,A,) (6) 
d t^oo 2at 

for example Bleher |5j uses this for d = 2. 

The dynamics may lead to anomalous diffusion, where the second moment may not be linear in t. 
Conventions differ as to how the diffusivity is generalised. Here we consider cases where the second moment 
is logarithmic. By analogy with the normal distribution case above we define 



lim i^i^ ^ x> . x> = ^tr2?,„-, ^ 7V(0, S„) (7) 

t^oo 2t\nt ^ d ^ tint ^ ' v ^ 



where the relevant limits exist, but now it is no longer clear whether Vij is related to S^-, as noted in }16[ . 
Convergence in distribution does not imply convergence of the moments, since the tails may decay slowly 
while spreading further with time; see Fig.[l] The authors of Ref. |16! numerically considered higher moments 
of A(n) /V n In n but it is likely that recent methods [35^ will lead to a proof of an anomalous exponent without 
a logarithm for these higher moments and with a logarithm for lower moments [34] ■ This proposed work 
does not include our case of interest (second moment) which is borderline between regimes of the normal 
and abnormal moments. We continue the discussion of this anomalous moment below. 

As with the normal diffusion, we assume that superdiffusion is possible in all directions; in general it is 
possible to have normal diffusion in some directions and superdiffusion in others, as discussed in 111) . 
Also, the definition for S makes sense with respect to discrete (collision) time n replacing t; we expect the 
relation S'^'^'^ = rS with t the mean free time; see 0111] and Sec.jgj Note however that 1?'''*'= is problematic 
since (A^) is infinite for a single free path as well as n free paths, and hence not susceptible to normalisation 
by any function of n. 

Noting Ai = /g Vi{s)ds we find (see Eq. 1.4 in Ref. 0) the expression 

l{\Aj)=tf {v,{0)v,is))ds- f s{v,{0)v^{s))ds (8) 
^ Jo Jo 

which leads to the usual continuous time Green-Kubo formula for Dij in the limit f — > oo, assuming the 
velocity autocorrelation function (VACF) decays sufficiently rapidly. If instead, we have {vi{0)vj{s)) ~ Cij/s, 
we find Vij = Cij . See Ref. [33] for rigorous results on continuous time decay of correlations for this system. 
Following Bleher [S], we note that in this (continuous time) case, the superdiffusion arises from the slow 
decay of the VACF, but that in the discrete case, the correlations decay rapidly, but superdiffusion arises 
from the infinite second moment of the displacement between collisions. 

We now return to the specific case of the periodic Lorentz gas and discuss existing results for = 2. The 
Lorentz gas is a billiard in an unbounded domain with a periodic collection of convex obstacles. If the free 
path length is bounded away from zero and infinity ("finite horizon"), the diffusion is normal, as noted in the 
Introduction. The Green-Kubo formula is valid, but there is no closed form expression for the correlation 
functions or diffusion coefficient. 
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Figure 1: The distribution of the scaled displacement A/ \/t In t for the square (d = 2) Lorentz gas with 
r = 0.4, showing convergence to the expected normal distribution, but with fat expanding tails at any finite 
time that support the argument in the text for an anomolously large second moment. Two features observed 
here at short times are oscillations due to the lattice, and spikes arising from ballistic orbits with no collisions 
that indicate the boundaries of the support of the distribution. 



In the infinite horizon case there are infinite strips that do not contain any scatterers, termed "corridors" 
that lead to logarithmic superdiffusion as described above. These corridors can be enumerated, leading 
to exact expressions for the superdiffusion coefficient. We assume that there are at least two non-parallel 
corridors, so that superdiffusion occurs in all directions. 

For a comparison of the coefficients existing in the literature, we need only consider the specific case of a 
Lorentz gas with a square lattice of unit spacing and disks of radius l/-\/8 < r < 1/2, for example the case 
r — 0.4 considered in numerical simulations below. This has only corridors in the horizontal and vertical 
directions. Bleher [5 following Eq. 8.5 writes a discrete superdiffusion coefficient (in our notation) 

,,2 



^'mchcr = £z (9) 



w 

where w = 1 — 2r is the width of the corridor, and 27rr is the circumference of a scatterer. As noted above 
the second moment of A is actually infinite. Bleher attempts to circumvent this issue by putting the n — > oo 
limit inside the average, but this also does not make sense at face value. If we make the interpretation 
pdisc ^ it agrees with the calculations presented here. 

Szasz and Varju ([IS], Eq. (2)), and later Chernov and Dolgopyat ([11 , Eq. (2.1)), give for the compo- 
nents of the covariance matrix 

^disc_Y- V'.(x)V,(x) (10) 



where cc is a fixed point of the collision map of the reduced (unit cell) space, assuming that each side of the 
corridor contains periodic copies of only one scatterer; there are four of these for each corridor, = l/(47rr) 
is the normalising factor for the invariant measure of this map, equal to half the reciprocal perimeter of the 
billiard, is the width of the corridor corresponding to x. Finally "^(x) is the vector parallel to the corridor 
giving the translation of x under the map in the full space, with components 'ipi{x). For the square lattice, 
ip{x) is a unit vector aligned with the coordinate axes, and taking account of the above mentioned factor of 
four, we find 

2 

"disc _ J^S- ■ 

Thus in contrast to the normal diffusion case it appears that 5 ^ 2T>, in fact E = D. This discrepancy is 
not noted in existing literature, presumably due to the different notations used. However, Balint, Chernov 
and Dolgopyat have reportedly proposed it for this (two dimensional) case, independently of and prior to 
the present author, and claim to have a proof [TU] . 

From Eq. ([s]) we can (heuristically) explain the discrepancy as follows. Given a random initial condition 
(according to the usual continuous time invariant probability measure /z), there is a probability ^ C/t that it 
will make a free flight for at least a long time t. This is much higher than the exponentially small probability 
of travelling this distance according to the normal distribution. Thus the latter can be accurate at best only 
to distances 

-a;2/(2Silnt) n 

Z,^^-— « - (12) 
\/2n^t\nt X 

that is. 



X « VEtln^t (13) 

Orbits dominated by a free flight much longer than this are well outside the converged part of the normal 
distribution and thus cannot contribute to its variance. Splitting the integrals in Eq. ([S]) into times less than, 
or greater than y/i, we find that the first term gives a contribution to D from the orbits that could contribute 
to the normal distribution of {C /2)tlnt, that is, exactly half its full value (modulo lower order terms like 
tlnlni arising from the logarithm in Eq. (13)). Thus we expect S = 2? in agreement with Eq. (11). Note 
that this argument depends on the specific mechanism for the slow decay of the VACF, and may not apply 
to all processes with logarithmic superdiffusion. It should however extend to the Lorentz gas in arbitrary 
dimension considered for the remainder of this paper. 

Thus we propose that the second moment is anomalous, specifically by a factor of two, giving the equation 
S = I? in two and higher dimensions. In d = 2 we agree with recent rigorous expression for S in Refs. |llll45j 
and Bleher's V. However, some of Bleher's discussion of P^isc g^^^^ ^ appears problematic. 
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Figure 2: Horizons in the cubic Lorentz gas, depicted with a 3 x 3 x 2 subset of the lattice, as viewed 
from one corner, for different scatterer radii, and ilhistrating the horizon terminology introduced in the text. 
If the scatterers are small enough that they don't touch (r < 0.5, left) there are planar horizons that are 
principal (and hence maximal) and not incipient. When the scatterers touch [r = 0.5, centre), they create 
an incipient principal horizon. When the spheres overlap (r > 0.5, right), cylindrical horizons are possible 
as seen in the figure; these are maximal but not principal. 

3 Summary of results 

As discussed previously for the case d = 2 O [HI |45] superdiffusion in an infinite horizon Lorentz gas is 
a result of the unbounded path length of the particle, thus an understanding requires close study of the 
properties of these paths. Here we extend the concept of a horizon (or "channel" ) as it appears in existing 
literature. Usual terminology defines free paths as bounded from above ("finite horizon") or infinite ("infinite 
horizon" ) with occasional variations ( "finite but unbounded" , also called "locally finite" , for some aperiodic 
Lorentz gases). Sanders j37j makes a finer classification of infinite horizon billiards in three dimensions using 
the dimensionality of the horizon (cylindrical or planar horizon billiards) . Here we define a horizon, not as a 
type of bound or dimension, but as a set of points whose trajectories have no collisions with the scatterers, 
and satisfy some natural conditions. 

Some definitions we need in this section (with more details in Sec. |4]) are as follows. We consider the 
compact configuration space C K'^/£ where the lattice C of covolume V is the discrete subgroup of 
defining the periodicity of the Lorentz gas. The scatterers cover a volume fraction V = 1 — \Mx\/V. The 
continuous time dynamics $* acts on the phase space M = x with My = S*^"^, the set of velocities of 
unit magnitude. The set Mt C M denotes initial conditions that do not collide with a scatterer for times in 
[0, t] . A horizon H C A/oo is of the form H = H^x Hy with Hy — §^'^^ n Vh and Vh C M'' is a linear space of 
dimension dn G {1 ... d — 1}, the "dimension of the horizon." The set is connected and invariant under 
translations in Vh ■ Projecting onto the d — dn dimensional orthogonal space Vj^ we obtain a lattice Cjj 
of covolume Vj^ and a set characterising the horizon H^. There is a natural invariant probability measure 
on M, denoted /i. The probability of surviving without a collision for time t is F{t) = fi{AIt), while the 
probability of remaining within a horizon H (or more precisely its spatial projection) is 

^HW-/^({(x,v):$^(x,v)ei?, xAf„, Vse[0,t]}). (14) 

We make further distinctions, illustrated in Fig. [2] A maximal horizon is one of highest dimension for 
the given billiard, a principal horizon is one of the highest dimension possible, d—1 (planar in Sanders' case) 
and an incipient horizon is one where H^^r has zero d — dn dimensional Lebesgue measure, hence Fnit) ~ 
for all t. For example, a principal horizon has a connected one-dimensional set that is an interval, of 
length wh, the "width" of the horizon. When this interval shrinks to a point, as by increasing the size of the 
scatterers, the horizon has scatterers touching it on both sides and becomes incipient. The set of maximal 
non- incipient horizons is denoted T-L. Lemma [T] in Sec. [5] asserts that this set is finite. 

We will see that the decay exponent of F{t) is determined by the dimensionality of the maximal horizon 
with special conditions related to incipient horizons, and 1?^^ and 'Df^^'^ are nonzero only for non-incipient 
principal horizons. 
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Sd 2 27r An 27r- -3- r((d+i)/2) 



^ 2 n ^ 2 r((d+2)/2) 

j_ 1 1 r{d/2) 

27r 4 2v/^r((d-l)/2) 



47r 



Table 1: Lebesgue measure of d-dimensional spheres, Sd and balls Vd, and the combination Gd = Sd-2/2S, 



d-l 



appearing in Eq. 27 below, including specific values of interest and the general formulae (last column). 



Although much progress has been made on higher dimensional billiards in the last decade, some major 
challenges remain ^44j . Existing results are not sufficient to provide a fully unconditional foundation for the 
study of diffusion in higher dimensional Lorentz gases, thus our results depend on the following conjectures. 

First, we assert that the leading order collision-free behaviour is determined only by maximal horizons, 
except where these are all incipient. This means in particular that incipient horizons, non-maximal horizons 
and overlaps between maximal horizons are all assumed to give subleading contributions. 

Conjecture 1. Consider a periodic Lorentz gas with at least one non-incipient maximal horizon. Then we 
have 

Fnit) ^ Fit) (15) 

Hen 

as t — ^ 00. 

The initial sentence of the conjecture implies that F{t) > for all time as required by the right hand 
side of the equation. The definition of periodic Lorentz gas considered in this paper is given at the start of 
Sec.H 

Remark 1. This conjecture describes only the horizons and not reflections from the scatterers. It may not 
be necessary for the scatterers to be dispersing, as its content is purely geometrical, and not dynamical. 
The main obstruction to proving it is probably the effect of maximal incipient horizons (see the discussion 
preceding Conj . [2] below) . 

Remark 2. Where there are only incipient horizons, it is likely that F{t) must decay to zero at finite time. 
We do not discuss this case here. 

Now we are ready to state the first main result: 

Theorem 1. Consider a periodic Lorentz gas satisfying Conjecture^ Then the free flight function F{t) 
satisfies 

Sd^,,J,^J,^A-%.,y)d.dy 

as t ^ 00. 

Here, A^'*'(x, y) is the "visibility" counting function, giving the (non-negative integer) number of ways x 
and y can be connected by a straight line entirely in i?^; for the specific examples considered here, this is 
either zero or one. Sd is defined in Tab. [l] 

Proof. Conjecture [T] asserts that only horizons in H need be considered to determine the asymptotic form 
of F{t). Lemma [1] (below) shows that the sum over H is finite, and lemma [2] (below) gives the value of each 
term in the sum. □ 

Sec. |6] gives the main application of this theorem, the expression ( [27| for F{t) in the case of the cubic 
Lorentz gas with principal horizons. 

If all maximal horizons are incipient, we cannot apply Conj.[l] However heuristic arguments and numerical 
simulations given in Sec. [7] suggest in the case of a principal incipient horizon, contributions from the infinite 
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number of horizons of next lower dimension can be summed, giving a t^"^ decay of F{t), as long as d < 6. 
For d > 6 the sum of contributions diverges, and it appears that F{t) has a decay intermediate between t^'^ 
and . Assuming a logarithm at the critical dimension d = 6 we can rather tentatively suggest: 

Conjecture 2. The free flight survival probability for a Lorentz gas with incipient (but no actual) principal 
horizon satisfies 

F{t) >i I t-^lnt d^6 (17) 

[ t-"" 1 <ad<2, d>6 

Finally we turn to the (super-)diffusion coefficient. Here, a further conjecture is required to establish 
estimates on velocity autocorrelation functions. Consider a flight of length t much larger than the lattice 
spacing; this has an angle from the horizon of order t~^. The particle makes a grazing collision with a 
scatterer and reflects with a typical angle of order t^^/^ so that the subsequent flight is typically of order 
^/i and further flights shorter still. Thus after any time proportional to t, we expect that many collisions 
have completely randomised the trajectories position relative to the scatterer and velocity direction, but not 
position in the full lattice. We state this as 

Conjecture 3. Consider a periodic Lorentz gas with flow $* and free flight function F{t) > for all t > Q. 
Let /, g : Af — > M denote zero-mean (with respect to /i) Holder functions. Then we have 

/ {f)igo^')dfi^o{F{t)) 

JM\Mt 

as i — 7> oo. 

Note that the \/t motivating the above conjecture holds for the Lorentz gas, but not for other models 
such as extended stadium billiards. Also, a connection between free flights and correlations for the three 
dimensional cubic Lorentz gas was proposed as early as 1984 19 . We can now address the diffusion properties 
of the Lorentz gas: 

Theorem 2. Consider a periodic Lorentz gas with non-incipient principal horizon. Assume Conjectures [7] 
ana Then the super- diffusion matrix is given by 

where n(_ff) is either of the two unit vectors normal to all velocities in H . 

Note that in contrast to two dimensions, it is preferable to describe horizons by their perpendicular rather 
than parallel directions, hence the projection operator appearing in this formula. 

Combining both theorems and taking the trace of the diffusion matrix, we find the simple relation: 

dV = lini tF{t) (19) 

independent of the details of the principal horizons. 

In the case that the principal horizons do not span the full space, the superdiffusion matrix is degenerate, 
and we expect that as in two dimensions, there will be normal diffusion in the direction(s) orthogonal to all 
these horizons if such motion is not completely blocked. Sec. [8] contains a numerical test of Conj.[3]and the 
proof of Thm. [2] with application to the discrete time diffusion coefficient given in Sec. [9] 

For both the free flight function and the diffusion coefficient, the limit of small r — >■ spherical scatterers 
leads to an increasing number of horizons. Taking a Mellin transform, both results can be formally related 
to the Riemann hypothesis (RH), the celebrated conjecture that the complex zeros of the Riemann zeta 
function C(s) lie on the line with real part 1/2 [T3]. The formal argument assumes that the contour integral 
can be closed avoiding the poles of C(s)~^ in the integrand in a convergent sequence of increasing semicircles. 
A complete proof of an equivalence with RH is not yet available even for the simpler related case of the open 
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circular billiard of Ref . \^ ; work on the latter is ongoing. Connections between lattice sums and RH is of 
course not new; see for example Here we make a connection between RH and transport in a definite 
physically motivated model; other physical systems with connections to RH are discussed in Ref. [41]. 

Further development of this subject could include verification and/or modification of the above conjectures 
on the theoretical side; the incipient case is particularly challenging due to the similarity with the example 
in Ref. |1] , hence almost certainly exponential growth of complexity. The case of multiple rationally located 
scatterers per unit cell naturally leads to Dirichlet twisted Epstein zeta functions, which may be interesting 
from a number theoretical point of view (see the discussion at the end of Sec. [6| . There are generalisations 
of the Lorentz gas which are almost periodic [HI [36] or have external fields [nftEl HZ] ■ 

On the applications side, periodic boundary conditions are commonly used in molecular dynamics sim- 
ulations, so a natural extension would be to investigate the (possibly unphysical) effects of these boundary 
conditions in simulations of low density systems of many particles. Transport in periodic structures with 
long free flights occur naturally in crystalline materials, nanotubes and periodically fabricated devices such 
as antidot arrays and various classes of metamaterials. For these, it should be noted that the soft potentials 
may modify the dynamics of particles moving close to a horizon direction. As with periodic time perturba- 
tions, the effect could be either more or less stable than free motion in the Lorentz gas depending on the 
details of the potential and the speed of the particle, however, as here, enumeration of horizons remains 
important for understanding transport properties. 

The remaining sections are structured as follows: Sec. |4] motivates and gives the necessary definitions 
related to the horizons. Sec. [5] develops the proof of Thm. [l] proving two lemmas and discussing Conj. [l] 
Sec. |6] then applies this theorem to principal horizons, giving F{t) ^ C/t with the constant determined 
exactly. Numerical methods are then discussed, and used to confirm these results. In the limit of small 
scatterers, connection is made with the Riemann Hypothesis. Sec.[7]then considers incipient horizons, where 
heuristic arguments and numerical simulations suggest Conj. [2] Sec. [8] then turns to the diffusion problem, 
proving Thm. [2] conditional on Conj s.[l] and [3| again for principal horizons. In this case the numerical results 
are inconclusive due to the logarithm being dominated by a constant at accessible numerical timescales. 
Finally, Sec. |9j treats the diffusion problem in discrete time. 

4 Definitions 

The phase space of the full periodic Lorentz gas is M = Mr^ x M„ where C M'' denotes the configuration 
space external to the scatterers and My = S''"^ is the set of velocities of unit magnitude. The same quantities 
without tildes are associated with the compact phase space of the dynamics considered modulo the lattice, 
so M = X My with = R'^/£\n. There are natural projections tt^ : M M^ and 7r„ : M -> M^. 
Here, £ is the lattice of periodicity of the Lorentz gas, ie a discrete subgroup of IR'' with finite covolume V, 
and in the main example given by Z''. Also, H, C M.'^/C denotes the union of the scatterers. The scatterers 
(connected components of $7) are assumed to be finite in number. The curvature of their boundaries has a 
positive lower bound, and the boundaries themselves are smooth. To allow for touching or overlapping 
spheres (for example), the scatterer boundaries may be non-smooth on a finite number of submanifolds of 
dimension at most d — 2; we require however that each scatterer be the union of a finite number of convex 
sets in this case. 

The uniform probability measure on M is = m/[V 5*^-1(1 — "P)] where m is the usual Lebesgue measure 
for M.'^ X S''"^, Sd-i is the Lebesgue measure of the sphere S''"^ (see Tab. [l]) and V is the packing fraction 
of the obstacles, equal to ma;(r2)/V where nix is Lebesgue measure on M.'^. 

The usual billiard flow is denoted : M — > M. Thus we have the free flight formula <I>*(x, v) = 
(x + vt mod £, v) in the absence of a collision. The time to the next collision is denoted r(x, v) = inf {t > 
: TTa; o $*(x, v) e dft}. If there is no future collision, we allow T(x, v) to be infinite. The set surviving for 
time t without a collision is Mt — {(x, v) G M : T(x, v) > t}, in terms of which F{t) = /i(Aft). One of our 
main aims is to characterise the large t asymptotics of this function. 

The flow 4>* is discontinuous at a collision, using the usual mirror reflection law, but we do not need the 
explicit formula here. Nor do we need to specify conventions for the flow at the exact moment of collision or 
for orbits that are tangent to the boundary or reach a non-smooth point on it; these sets are of zero measure. 
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In three or more dimensions, collections of orbits with large or infinite T are of different types. We now 
formalise the terms introduced in the previous section. We want to characterise subsets of Moo , ie collections 
of points such that T = oo. 

The free dynamics is uniform motion on a torus. Thus a trajectory starting from a point in M must 
either be dense (hence have finite T) or move parallel to a linear subspace spanned by vectors in C. Thus 
for each x e tTxMoo there are one or more (possibly infinitely many) linear subspaces of directions in which 
there are no collisions. Let us choose one, Vh which is not contained in another such subspace. Then, we 
can define 

= y^ns^^i (20) 

Hx = {x e : (x, H,) G Afoo} (21) 

The equation defining free motion ensures that Hr,. is invariant under translations in Vh- Thus we can 
decompose it as 

Hx = H^(S Vh/Ch (22) 

where Ch —CD Vh is a sublatticc of the original lattice C. Initial conditions contained in and having 
velocities close to vectors in Hy will survive for a long time before colliding, as long as they remain in a 
connected component of H^- This motivates the following definition of a horizon: 

Definition 1. A Horizon is a set _ff C A/qo which satisfies 

(a) H is sl direct product of sets Hr^ C and Hy C My, 

(b) H is inextensible, ie not a proper subset of another horizon, 

(c) Hx is connected. 



The above discussion implies that all horizons are of the form given in Eqs. (20-22) above, with the 
additional condition that Hj. be a connected component of H^., that is, 

H = HxX Hy (23) 
H.. = H^® Vh/Ch (24) 
Hy = VhC]^'^-^ (25) 

The dimension of the horizon dn e{l...(i — l}is defined to be the dimension of V/f, which is also the 
dimension of Ch- When dn = 1, Hy consists of two opposite points, so H itself is not connected, even 
though Hx is connected. When dn > 1, is also connected. 

An example where H^ may have more than one connected component is that of a square lattice, ie 
£ = Z^, with a circular scatterer of radius ri = 0.2 at the origin and an additional scatterer of radius 
r2 = 0.1 at (1/2, 1/2). Particles moving horizontally will not hit any scatterer if their y coordinate lies in 
either of the intervals (0.2, 0.4), (0.6, 0.8). Thus there are two horizons with the same Vh in this case. 

We denote the d — dn dimensional covolume of Ch in Vh by Vh- The space of dimension d — dn 
perpendicular to Vh is denoted Vh - The original lattice induces a lattice Ch C V^, given by elements of 
C/ Ch projected onto Vh along Vh- This has covolume = V /Vh since a unit cell of the original lattice 
can be formed by the product of two unit cells of lower dimension in orthogonal subspaces. In the original 
space M each horizon H corresponds to a periodic collection of connected components; lattice translations 
in Ch preserve each component, while lattice translations in C\C-}i map one component to another. 

When d = 2 we can have only horizons with dimension dn = 1, corresponding to the channels in previous 
literature (see Sec. [2]). When d = 3 we can have a "planar horizon", dn = 2, in which Hy is a great 
circle in S^, a "cylindrical horizon", dn = 1, in which Hy consist of two symmetrically placed points; this 
terminology comes from |37| . More than one type of horizon may exist in the same billiard. For example 
if the scatterers do not cross the boundaries of the unit cell in three dimensions, there are planar horizons 
along the boundaries; there may also be a horizon of cylindrical type in the interior. Finite horizon, in 
which no horizons exist, are also possible in any dimension d, and becomes increasingly probable for random 
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arrangements of scatterers of fixed density in a unit cell of increasing size. Nevertheless no explicit example 
of a finite horizon periodic Lorentz gas in d > 2 has yet been shown. 

Recall the discussion of the free flight function, F{t) above. We can similarly define the probability of 
remaining for at least a given time in the spatial projection of a horizon, F}i(t)\ this definition has already 



been given in Eq. ( 14 ) . It is useful to make some definitions related to the dimension: 



Definition 2. Maximal horizon: A horizon of dimension greater than or equal to all other horizons of the 
given billiard. 

Definition 3. Principal horizon: A horizon H of the highest possible dimension, that is, dn — d — \. 

Clearly a principal horizon is also maximal. Anomalous diffusion properties arising from a. 1/t decay of 
velocity autocorrelation functions are expected only for principal horizons; see Ref. [331 137] and below. 

We now make an important observation: For dimensions d > 2 we have a qualitatively new type of 
horizon, corresponding to the border between horizons of one dimension and the next. 

Definition 4. Incipient horizon: A horizon for which has zero d — du dimensional Lebesgue measure. 
An incipient maximal or incipient principal horizon will denote a horizon that is incipient and also maximal 
or principal, respectively. 

When d = 2, an incipient horizon is an orbit tangent to scatterers on both sides, hence isolated and 
having no effect on the dynamics; this is a rather trivial form of horizon. In dimensions d > 2, a incipient 
horizon occurs when for a horizon of dimension dn > 1 shrinks to a zero measure set, so that Ffjit) = 0. 
This is easiest to visualise in the three dimensional cubic Lorentz gas, when the planar horizons shrink to 
zero width touching scatterers on both sides as the radius approaches 1/2 from below; see Fig. [2] 

Within an incipient horizon H ^ there is a lattice of points Cine in where the scatterers touch the 
horizon. The set of velocities within Hy that do not come arbitrarily close to the scatterers consists of those 
contained in horizons for a Lorentz gas of dimension du with infinitesimal scatterers placed at points of Cine- 
Each of these velocities corresponds to a horizon of dimension < dn — 1 in iJ and also in the original space 
M. However, unlike in the Lorentz gas with finite sized scatterers, the number of horizons of the highest 
dimension dn — 1 is now countably infinite, parallel to all possible du — 1 dimensional subspaces spanned by 
vectors in Cine- These become finite in number as soon as the radius of the scatterers is increased and the 
incipient horizon destroyed. Thus an incipient horizon of dimension dn > 1 is in general associated with an 
infinite number of horizons of dimension dfj — 1. This leads to the interesting effects discussed in Conj. [2] 
above and Sec. below. 

We also remark that incipient horizons contain orbits tangent to infinitely many scatterers from more 
than one direction, like the very recent example of exponential growth of complexity found in finite horizon 
multidimensional billiards [4]. Here too, the singularity structure in the vicinity of such orbits is extremely 
complicated. One final definition: 

Definition 5. % is the set of non-incipient maximal horizons. 



5 Two lemmas 

We now consider the first two propositions in section 2 of Bleher [5] in the more general case of arbitrary 
dimension, which are necessary for the proof of Thm. [l] The first proposition can be naturally generalised, 
and is given as the first lemma here. The second proposition is not valid in d > 2, and must be replaced by 
Conj. [T] above. The contribution from each horizon, also needed for Thm. [l] forms the content of the second 
lemma. 

We generalise Bleher's [5] Prop. 2.1: 

Lemma 1. For any periodic Lorentz gas with infinite horizon there are finitely many maximal (possibly 
incipient) horizons. 
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Proof. Assume there are an infinite number of maximal Irorizons, eacli with its velocity subspace Vh- The 
set of Vh is either finite or infinite. If it is finite, at least one Vh corresponds to infinitely many H, so that 
Hx has infinitely many connected components. But this is not possible since the scatterers are comprised 
of only a finite number of convex pieces. Thus we conclude there are also an infinite number of Vh, leading 
to unbounded d//-dimensional covolumes Vh- The horizons thus fill at least one space of dimension greater 
than dn more and more densely, which is only possible if the intersection of the interior of the scatterers 
with this space is empty. Thus the space forms its own (possibly incipient) horizon with a dimension greater 
than dn, violating the assumption that the original horizons are maximal. □ 

This lemma ensures that the sum over maximal horizons appearing in Thm.[l]is finite. 

Bleher's p! Prop. 2.2 asserts that in any sufficiently long free path (greater than a fixed constant Rq), 
all but extreme parts of bounded (less than another constant total length lie in the space projection 
of a horizon H, or in the terminology of that paper, a corridor. For d > 2 this is not true, since there is in 
general a complicated set of horizons of differing dimension connected in a nontrivial manner. 

All we can hope for is conjecture [l] given in the introduction, which we now discuss in more detail. The 
sum over H £ H ignores lower dimensional horizons, some of which are connected to the maximal horizons. 
It also double counts overlaps of the maximal horizons. For example, in the three dimensional Lorentz gas, 
with a cubic arrangement of spherical scatterers of radius r < 1/2 (considered in more detail later), there 
are planar horizons parallel to the {x,y), {x,z) and (y, z) planes, for example the set {{x,y,z) G M'^ : r < 
z < 1 — r}. Any two of these intersect in cylindrical sets parallel to one of the coordinate axes. There are 
also lower dimensional (ie cylindrical) horizons that are not part of the planar horizons, for example the set 
{{x, y,z) G : < x,y < r,x^ + > r^}. Conjecture [l] asserts that these lower dimensional effects do not 
affect the leading order free flight function F{t) if there is a non-incipient maximal horizon. 

The exclusion of purely incipient maximal horizons is necessary, as in this case Fnit) — 0. For dn > 1 
this is associated with an infinite number of horizons of lower dimension, whose contributions (including 
intersections) would need to be carefully enumerated to calculate F(t)] we give some heuristic arguments 
and numerical simulations for the case of incipient principal horizons in Sec. [Tj 

We now come to the most important result of this section: 

Lemma 2. For a non-incipient horizon H , we have as t oo 

S,^^,J J A--i.,y)d.dy 

Proof. We split v into components in and Vh subspaces: v = + v", and similarly for x. The 
trajectories not reaching the boundary of the horizon before time t are given by G (if^"^(x) — x)/t to 
leading order in t. Here, H^^^{x.) is the set of points "visible" from x in ie the set connected by a 
straight line that does not leave H;^. The tilde indicates that this set is defined on the full Vj^, not modulo 
the lattice, ie more than one lattice translate of a point may be visible to a given x G Hx- 



The definition of Fnit) is given in Eq. (14 1. Now is just given by m/[VSd-i{l — V)] where V is the 
packing fraction of scatterers and m is Lcbesgue measure. We integrate over all x" which is the covolume of 
Ch in Vh, Vh- The measure over x-*- and v-'-, multiplied by , is given by the integral 



which can be written more symmetrically as 



iijr(x)|dx 



A™(x,y)dxdy 



Since is of order t^^, v" remains of unit magnitude to leading order, and its measure is simply Sd„-i- 
Taking account of the normalisation factors and the relation V = VnVjj we get the result stated in the 
lemma. □ 
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6 Principal horizons and the free flight function 

In the case of a principal horizon, we note that is one dimensional, so that the double integral in Lemma 
2 is simply Otherwise i?^ is typically non-convex and the integral is not easy to evaluate (an example 

is provided in the next section). 

Possible principal horizons can be enumerated using the dual lattice £* = {u e M'^ : u-v e Z'', Vv e £} 
where u-v denotes the usual inner product in M.'^. In particular, for each lattice Ch of dimension d—1, either 
of the two vectors perpendicular to Vh and of length (V^)~^ lie in £* . It is not hard to see that there is 
a 1:1 correspondence between codimension-1 lattices and primitive elements of C* modulo inversion, where 
primitive indicates that it is not a multiple of an element of C* by an integer other than ±1, and inversion 
is multiplication by —1. Note that the zero vector is not primitive. We denote such a vector by 1 = n/^^/V^ 
with magnitude L = |1| = (V^)~^ so that is a unit vector perpendicular to Vh- The subscript H is 
suppressed for the vectors 1 and their magnitudes L. For each principal horizon there is an ambiguity in the 
sign of 1 and uh that does not affect any of the results. When summing over all primitive lattice vectors, 
the set denoted C* , we need to divide by 2 to take account of this inversion symmetry. 

Given a vector 1 S £*, a principal horizon will exist if the perpendicular projection of the scatterers onto 
the line parallel to 1 does not cover the whole of that line. If there is one spherical scatterer of radius r 
per unit cell, this corresponds to the condition that the width of the horizon wh — Vjj — 2r — L^^ — 2r is 
non-negative. Thus for a given r, there are a finite number of principal horizons corresponding to vectors in 
C* modulo inversion with lengths less than (2r)~^; this number grows without bound as r approaches zero. 
For a finite number of arbitrary convex scatterers, this argument gives an upper bound if 2r is interpreted 
as the smallest width of any of the scatterers. 

Returning to the case of a single spherical scatterer of radius r per unit cell, we have V = V^r'' (see 
Tab.[l]). Using Thm.[l| we find that 



where ^2(2;) — x'^u{x), u{x) the unit step function, equal to 1 for positive argument and zero otherwise, and 
L — |1| as above. The geometrical factor Gd = Sd-2/ {2,Sd~i) is given in Tab. 1, where the 2 comes from the 
inversion symmetry. 

Numerical simulations were performed using the author's C-|— I- code, which is designed to simulate 
arbitrary piecewise quadric billiards in arbitrary dimensions. The flight time is thus given the solution of 
a quadratic equation; simplifications for the cases where the matrix defining the quadratic form is zero or 
the unit matrix are used for greater efficiency. The main numerical difficulty with billiards is to ensure that 
the particle remains on the correct side of the boundary, while retaining as much accuracy as possible for 
almost tangential collisions. These can be addressed by disallowing very rapid multiple collisions within a 
very small perpendicular distance of the boundary. Tests were made comparing 16 and 19 digit precision 
calculations, and by an explicit investigation of almost tangential collisions. 

For the calculations presented here, the spherical scatterer is contained in a cubic cell. On reaching the 
boundary, the particle is transferred to the opposite side, updating the local position while keeping track 
of the overall displacement. The main obstacle to simulations in dimensions above d — 10 considered here 
is that the volume fraction of balls decreases rapidly with dimension for the cubic lattice. Other lattices 
(for example the close packed £'§ lattice) may be considered in the future, however a higher volume fraction 
comes at the cost of a more complicated unit cell. It is difficult to get an intuitive grasp of the structure 
of horizons of non-cubic lattices above three dimensions; one approach would be to numerically infer them 
from F{t) via Thm. [l] 

Numerical simulations of F{t) were performed for the case of a cubic lattice with r — 0.4 where there is 
only a single contribution in the above sum, for 2 < d < 10 in Fig. [3] Agreement is good for low dimensions, 
while for higher dimensions it is plausible but limited by finite statistics together with longer times required 
for convergence to the final C /t form. In the case of r = 0.6 the scatterers are overlapping, and the maximal 
horizons have dimension d~2. Here we expect C/i^ behaviour for d > 3 but do not have an explicit analytic 
formula; see Fig. [4] which shows good qualitative agreement for small dimensions, and slow convergence for 




(27) 
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Figure 3: The free flight function for dimensions 2 < d < 10 for r = 0.4 showing C /t behaviour. The inset 



shows comparison with the predictions of Eq. (27); here the numerical results are dotted and the predictions 
the solid straight lines. 

higher dimensions. The intermediate case of r = 0.5 (incipient principal horizon) is considered in section [t] 
below. 

We now seek the small r asymptotics of the above expression. First we remove the primitivity condition 
using Mobius inversion, 

^ CO ' 

m - ,n - L^) ^ ^ ^.{C)CLu,{{CL)-^ - 2r) (28) 
' c=i lec 

where the prime indicates that the sum excludes the zero vector. We take out the geometric factor and use 
a Mellin transform on the rest to extract the small r behaviour 

p.,^ Gg 1 Ar-^E{{s + l)/2;C*)ds 

where E{s'; C*) is the Epstein zeta function for the dual lattice, 

The purpose of the Mellin transform is to express the contour integral as a sum of residues by Cauchy's 
theorem, closing the contour to the left. This sum over poles directly gives an expansion in powers of r 
(also logarithms if the poles are not single). It naturally depends on being able to show that there exists 
a sequence of roughly semicircular contours for negative real part of s and avoiding poles of the integrand, 
with integral converging to zero; we are unable to prove this here. 

The integrand includes the Riemann zeta function C(s) in the denominator, so the powers of r appearing 
in the residue sum depend on the locations of the zeros of C(s). The celebrated Riemann Hypothesis [1^ 
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Figure 4: The free flight function for dimensions 2 < d < 10 for r = 0.6 (overlapping) scatterers, showing a 
collision at bounded time for d = 2 and approach to C/f^ for d > 2; the rate of convergence with t is slow 
for the higher dimensions considered. 



asserts that the complex zeros lie on the line with real part of s equal to 1/2; the zeros on the real axis are 
all known. 

We also need some analytic properties of the Epstein zeta function in the numerator; see |38j . The function 
E(z,C*) for arbitrary lattice of full dimension C* has an analytic continuation in the whole complex plane 
except for a single pole a,t z = d/2 with residue Sd-i/2, a special value E{0; £*) = —! and zeros at all negative 
integers. The following statements [13] show how it generalises the Riemann zeta function C(s): Considering 
the (self-dual) cubic lattice C = C* = Z'', we have E{s';Z) = 2C(2s') and E{s';I?) = ^C,{s')l3{s') where /3 is 
the Dirichlet /3-function, that is, the nontrivial Dirichlet i-function modulo 4. In four dimensions we have 
i?(s';Z^) — 8(1 — 2^^^'' )C(s')C(s' — 1), which is already enough to show that the Epstein zeta functions do 
not all satisfy a Riemann hypothesis with all complex zeros on a single vertical line in the complex plane; 
note that Ref. [43J has 2^~* (a misprint). There is no factorisation for the physically interesting case of 
three dimensions, as for example there are lattice vectors with lengths a/S and a/S but not the product vTS. 

The above integral has poles ats = d— 1,— 1,— 2 and at complex values with real part —1/2 due to the 
C(s + 1) in the denominator (also elsewhere if the Riemann Hypothesis is false). Assuming that the contour 
can be continued to the left, avoiding poles of the Riemann zeta function, we find for the d-dimensional cubic 
lattice, 

= 2^dr((/+3V2K(d).^-^ + 

where the leading coefficient is the same as found for the Boltzmann-Grad limit in |31j . and the next term 
is valid for all (5 > if the Riemann Hypothesis is true. Numerically this term is small in magnitude for r 
not too small; the next term 

Sd-2Vd of< 

2d'i{d^ ^d)C{d) \ 

is quite visible; see Fig.[5j This argument with the contour is similar to that of the survival probability of the 
circular billiard, Ref. [5], and as in that case suggests an equivalence of Eq. 31 to the Riemann Hypothesis, 
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Figure 5: The sum in Eq.(27) minus the leading term in Eq. (31), compared with the series obtained by 



summing the residues of poles on the real axis in Eq. (29), the smooth curve. The smooth terms give a good 



approximation for the range of r considered; for smaller r the complex poles coming from the zeta function 
dominate. 



requiring a further careful argument. 

In the case of two or more scatterers per unit cell with rationally located centres and common denominator 
q, the same calculation can be performed. The Mellin transform then gives an expression of the form 



Ls+i 



E 



Cs+l 



(32) 



where x is a Dirichlet character mod q, f^is) and m are explicitly known, the sum over C can be expressed 
in terms of known Dirichlet L-functions, relating transport in this Lorentz gas to the generalised Riemann 
Hypothesis. The sum over 1 gives Dirichlet-twisted Epstein zeta-functions. The latter may be difficult to 
characterise in d = 3 due to the lack of a factorisation of the untwisted Epstein zeta function. 



7 Incipient horizons and the critical dimension 

The results of the last section do not apply directly to the case of incipient principal horizons, as their width 
is zero, and they are associated with an infinite number of horizons of the next lower dimension. There is 
however at least one useful piece of information to be gained, which we now describe. 

Lemma [2] shows that the contribution to F{t) from a horizon of dimension d — 2 is of order t^^. A finite 
number of these contributions are thus of the same order, however a more careful calculation is required for 
the present case, where there are an infinite number of contributions. 

We consider a single example, but one which is likely to have more general implications. Consider the 
Lorentz gas with a cubic lattice of arbitrary dimension d > 2 and r — 1/2. The scatterers just touch, leading 
to incipient horizons in the coordinate planes and no other principal horizons. In one of these planes, for 
example at — 1/2 for the d-th coordinate, the cross-section of the scatterers are reduced to the points of 
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Figure 6: The region defined by Eq. [33] 



2<i-i j-Q f^i^Q other d—1 coordinates, equivalent to the problem in the previous section, with dimension d—1 
and r = 0. Thus there are horizons of dimension d—2 corresponding to all non-zero lattice vectors 1 G Z'^^^ 
modulo inversion. For any given lattice vector, there will also be a finite extent in the d-th coordinate 
direction, which we now quantify. 

The space is two dimensional, in the plane spanned by 1 and e^, the unit vector in the d-th coordinate. 
Projected orthogonally into this plane, the lattice has basis vectors (hence of magnitude L^^, where 
L = |1|, consistent with previous notation) and e^. This means consists of the area confined between four 
circles of radius r — 1/2 and centred at (0, 0), {L^^, 0), (0, 1), {L^^, 1) with respect to the orthonormal basis 
{l/L,ed}; we will denote coordinates with respect to this basis as {x,y). Pairs of scatterers separated by 
L^^ overlap (except when L = 1) and pairs displaced by 1 just touch. The region not within radius r — 1/2 
of the centre of any scatterer is contained in the rectangle defined by < a: < L^^ and — L^^ < U < 
1 — \y/l — In the limit as L — > oo the height in the y direction is L~^/2 to leading order. In this limit, 
we scale in both directions, x — L~^/2 ~ L^^S, and y— 1/2 = L~^ri and keep only the highest order nontrivial 
terms in the equation, and find that is defined by a limiting shape bounded by four parabolas. 



\v\<m-m' 



ICI < 1/2 



(33) 



This region is plotted in Fig. [6] There is thus a limiting nonzero probability for two points in the region to 
be mutually "visible" , ie for which the line joining the points remains in the region. 

The double integral over in Lemma 2 is thus asymptotic to a constant (numerically found to be 

approximately 0.02746) multiplied by L^^ as i — > oo. The extra factor of L in Lemma^lshows that 



FH{t) 



(34) 



for some constant Cd containing only geometrical factors, in the appropriate limit (first t — ^ cx), then L — >■ oo). 
Summing over 1 S gives the Epstein zeta function i?(5/2; Z''~^)/C(5). This is a convergent series for 

d < 5 and divergent for d > 6. Thus d = 6 is a critical dimension for these incipient horizons. The divergence 
of the sum indicates that F{t) has a slower decay with t than t^^ but a detailed calculation is required. 
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Figure 7: The free flight survival probabihty for dimensions 2 < c? < 10 for r = 0.5 (touching) scatterers, 
showing escape at finite time for d — 2, approach to C/t'^ for 2 < d < 5 and slower decay for d > 6 as 
predicted by theory; fits for the relevant exponents are given in Tab. [2j 



d 


3 


4 


5 


6 


7 


8 


9 10 


n 


= 10« 


2.07 


2.17 


2.19 


1.99 


1.82 


1.70 


1.64 1.58 


n 


= 109 


2.02 


2.20 


2.17 


1.98 


1.80 


1.67 




n 


= 10^" 


2.03 


2.08 


2.14 


1.95 









Table 2: Numerical estimate of the exponent a from Conjecture |2] for cubic Lorentz gases with r = 1/2, 
sample size n, and fitted for free flight survival probabihties 10^/n < F{t) < 10''/n. 

taking into account the overlap of the relevant horizons, to yield a precise form. It is unlikely that the result 
would reach that of the non-incipient principal horizon, that is, t~^; we propose Conj. [2] above. 

Numerical tests of this conjecture are somewhat difficult due to the low packing fraction in high dimensions 
(so that individual trajectories survive for long times) and slow convergence to the long time limit. Scatterer 
radii r — 0.4 and r — 0.6 were discussed previously; see Figs. |3] and |4] For a scatterer of radius 0.5 and 
dimension d > 6 the exponent itself is unknown; results are shown in Fig. [7] with fitted exponent in Tab. 2, 
and show exponents greater than 2 for low dimensions due to the finite times involved. 

8 Correlations and diffusion 

The above expression for F{t) is also relevant to correlation functions, since conjecture[3]asserts that following 
a long fiight, correlations decay very rapidly. For describing diffusion, we need velocity correlation functions. 
Figs, [s] and |9] show the velocity autocorrelation function (used below) compared with the sum of horizons 
predicted by Conj. |3]in the case of a three dimensional cubic Lorentz gas. 

Recalling the discussion in previous sections, we are interested in the case of nonincipient principal hori- 
zons for which the function F{t) decays as C/t; only in this case do we expect Eq. ([8| to give superdiffusion, 
ie a nonzero value of V. As discussed for F{t), however, the coefficient is explicit. Now we give the proof of 
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Figure 8: For r < 0.5, Eq. 27 gives the limit of tF{t) as t — > cx) as a sum over horizons, assuming Conj. [T] 
Conj. [3] apphed to the components of the velocity asserts that the same limit is given by f(vo • Vj) related 
to the diffusion by Eq. ([s]). Here this correlation function is shown with a fit to the form Cr/t, for the three 
dimensional cubic Lorentz gas. 



100 




0.01 



Figure 9: See Fig. [Sj Comp arison of the fitted Cr (single points) with limt^aa tF{t) from Eq. (smooth 
curve) . 
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Thm. [2] 



Proof. Assuming conjecture [3j the velocity correlations are determined to leading order by the integral over 
orbits that do not reach a scatterer in the given time, that is, Alt- For these orbits, v is a constant, so 
conjecture [T] can be applied. We then use the same method as in the proof of Lemma 2, except that the 
integrand varies significantly over the v" space, ie Hy. Given that the perpendicular space is one-dimensional, 
we find a simple result in terms of the width of the horizon wh defined at the beginning of Sec. [6] 

{vmvM - E ' ' (35) 

as t — !■ oo. Now Hy consists of all unit vectors v perpendicular to a unit vector n, which depends on 
the horizon H but we will suppress this in the notation. This n is parallel to a vector in C* as discussed 
previously. Switching to a coordinate system based on n we can evaluate the integral to give 

ViVj = {Sij - ninj)Vd-i (36) 

The cases d> 2 and d = 2 need to be derived separately, but give the same result. Putting this together we 
find the stated result. □ 

For the case of the cubic lattice, the superdiffusion matrix is isotropic. For l/(2\/2) < r < 1/2 we have 
only d horizons of the form (1, 0, 0, . . .) and find that 2?,%,- = VSij with 

(1 - 2rf rid/2) 

while for l/(2\/3) < r < 1/(2^2) there is an additional term coming from the d{d~^ 1) horizons of the form 
(1, ±1, 0, . . .), which is thus larger by a factor of order d. In general this means that the large d asymptotics 
at fixed r has a different form every time a new horizon is introduced. 

For the limit of small scatterers and fixed d the sum over horizons is exactly the same as for F{t) in 
and we can directly combine Eqs. ( |19|31[ ) to get 



^- 2^ci^r((/r3V2)c(^).^-^ +^(-''^") 

for (5 > 0. Again it appears that this statement is equivalent to the Riemann Hypothesis under the same 
caveat (convergence of a sequence of contour integrals) as in Sec. [6j 

We can also compute the diffusion numerically, however even for the lowest dimensions, normal diffusion 
unrelated to the horizons dominates at accessible timescales; see Fig. [lO| 



9 Diffusion in discrete time 

In order to compare our explicit expressions with other literature, we note that [IS] uses exclusively discrete 
time properties (using n rather than t in the definitions), while O link discrete to continuous time 
preperties. Recall the discussion in Sec. [2} 2?'^"^'= is formally infinite, however we can discuss the covariance 
of the limiting distribution, which is related to the superdiffusion coefficient by S = 2? and to its discrete 
time version by S'^"^'^ = tS. The mean free time r is given by [5] 

- (39) 



\^VLC^E\ Vd- 



Here \dVl n i?| is the measure of the boundary of the scatterers, equal to Sd-ir'^ ^ in the case of a single 
scatterer per unit cell, without overlapping. Thus in the latter case we find 

^ 1 - Frf r'^ 



(40) 
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Figure 10: The superdifFusion coefficient V numerically (dotted lines) and theoretically (horizontal lines). 
At accessible timescales normal diffusion processes dominate. 

which when combined with Thm.[2]and the above relations gives a simpler prefactor: 
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